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We compute the tree-level connected four-point function of the primordial curvature perturbation 
for a fairly general minimally coupled single field inflationary model, where the inflaton's Lagrangian 
is a general function of the scalar field and its first derivatives. This model includes K-inflation and 
DBI-inflation as particular cases. We show that, at the leading order in the slow-roll expansion and 
in the small sound speed limit, there are two important tree-level diagrams for the trispectrum. One 
is a diagram where a scalar mode is exchanged and the other is a diagram where the interaction 
occurs at a point, i.e. a contact interaction diagram. The scalar exchange contribution is comparable 
to the contact interaction contribution. For the DBI-inflation model, in the so-called equilateral 
configuration, the scalar exchange trispectrum is maximized when the angles between the four 
momentum vectors are equal and in this case the amplitude of the trispectrum from the scalar 
exchange is one order of magnitude higher than the contact interaction trispectrum. 



I. INTRODUCTION 
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Precise measurements of the cosmic microwave background (CMB) anisotropies such as those obtained by the 
WMAP satellite [l[ provide valuable information on the very early universe. Any theoretical model that attempts 
to explain the evolution of the universe before the big bang nucleosynthesis will also have to explain the observed 
CMB anisotropies. The power spectrum of these primordial anisotropies is nearly scale invariant and it contains 
almost all the information on primordial perturbations. However, a small amount of non-Gaussianity is still allowed 
by the data 0,0]. 

The information contained in this non-Gaussian component will contribute to a huge advance in 
our understanding of the very early universe. For example, the simplest and most popular inflationary model based 
on a single field with a canonical kinetic term satisfying slow-roll conditions and standard initial conditions predicts 
that the level of primordial non-Gaussianit is actually unobservably small [4|, even with future CMB experiments like 
the Planck satellite [5j . 

Recently however, there are some hints of non-Gaussianity 0, H, 0] in the CMB anisotropies. These are one (or 
so) standard deviation hints but they motivate tremendous efforts from the theoretical cosmology community to 



create inflationary models that can produce large levels of non-Gaussianity both seen in the bis 
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can produce sizeable non-Gaussianity. They basically relax one (or more) of the following conditions: 
slow-roll, standard kinetic energy and standard initial conditions. 

In the present work, we will consider a single field model with a non-canonical kinetic term. It has been shown 
[IH EH H3, HE] that this kind of models predicts large bispectrum and also large trispectrum of the primordial 
curvature perturbation. In particular the DBI-inflation model (54[, which is constructed within string theory, is such 
an example. For current and stringent observational constraints and consequences of DBI-inflation see [H, HH H3, HI, 

m miH, mm. 

We will focus our attention on the trispectrum of the primordial curvature perturbation. The current observational 
bound [64, 65] on the trispectrum is rather weak |tjv.l| < 10 8 , where tnl parameterizes the size of the trispectrum, 
but in the near future with the Planck satellite it will improve to \tnl\ ~ 560 [6(|. We should note that these bounds 
depend on the shape of the trispectrum [67j and the values quoted here can be applied only for local models of 
non-Gaussianity. Thus it is important to calculate the exact form of the trispectrum in the DBI-inflation model and 
construct an estimator to measure the trispectrum of the CMB anisotropies. 
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In [68(, Seery et al. calculated the contact interaction trispectrum for slow-roll multiple fields inflation but with 
standard kinetic terms. They showed that at horizon crossing the contact interaction trispectrum is too small to be 
observed. They found r^f ~ r, where r is the tensor to scalar ratio, that is currently constrained to the range r < 0.22 
at 95% confidence level Recently, Seery et al. [69| , computed the graviton exchange trispectrum and they showed 
that the above conclusion does not change, i.e. the total non-linearity parameter t^l coming from both the contact 
interaction and the graviton exchange trispectrum is still unobservably small and of order r. Also recently, Gao and 
Hu [49| computed the leading order trispectrum from entropy perturbations in multifield DBI-inflation model. 

In [25j (see also [30|), the authors computed the trispectrum at leading order in slow-roll and in the small sound 
speed limit for a fairly general inflationary model, where the inflaton's Lagrangian is a general function of the scalar 
field and its first derivatives. They only calculated the contact interaction trispectrum and they argue that for certain 
models (DBI-inflation inclusive) the trispectrum can reach values in the observable range of the Planck satellite. 

In the present work, we will calculate the scalar exchange trispectrum for the same model. We will show that the 
scalar exchange trispectrum is of the same importance as the contact interaction trispectrum. This completes the 
work of [25| (see also [3(3]) and gives the final answer that should be constrained when comparing the trispectrum for 
models with large non-Gaussianity with observations. 

The rest of this paper is organized as follows. In the next section, we shall introduce the model and some basic 
notations. In section HTI1 we will consider non-linear perturbations in the action up to quartic order. In subsection 
MI Al we will present a simple argument to prove that the amplitude of the scalar exchange trispectrum is of the 
same order of magnitude as the contact interaction trispectrum. In subsection llll B[ we will calculate the trispectrum 
of the primordial curvature perturbation at leading order in slow-roll and in the small sound speed limit when only 
contact interactions are taken into account. In subsection MI Cl we shall calculate the scalar exchange trispectrum. 
At the end of this subsection, we present for the first time the complete trispectrum at leading order in the previously 
mentioned approximations. Section IIVI is devoted to the study of the shape of the trispectrum in the so-called 
equilateral configuration for the DBI-inflation model. We will also calculate the non-linearity parameter t^l for 
different models. The conclusions of this work are described in section [Vj 



II. THE MODEL 



In this section, we will introduce the inflationary model under study. We shall present the background equations 
of motion and define the slow-variation parameters. At the end of the section, we show the results for the primordial 
power spectrum of the curvature perturbation and the spectral index [70| . 

We will consider a fairly general class of models described by the following action 

S = \ J d 4 x^—g [M 2 Pl R + 2P(X, 0)] , (1) 

where <j> is the inflaton field, Mpi is the reduced Planck mass, R is the Ricci scalar and X = ~(l/2)g^ d^dycf), where 
g^ v is the metric tensor. We label the inflaton's Lagrangian by P and we assume that it is a well behaved function of 
two variables, the inflaton field and X. Throughout this work, we use a system of units where the Planck constant H, 
the speed of light c and the reduced Planck mass Mpi are set to unity. This general Lagrangian includes as particular 
cases the DBI-inflation model [IHHil and the K-inflation model [7l|. In the DBI-inflation model P(X, <ft) is given by 



P(X, cj>) = -f{ct>)- Vl - 2/(^)1 - (2) 

where /(</>) and V(4>) are functions of the scalar field determined by string theory configurations. At the background 
level, we assume that our Universe is well described by a flat, homogeneous and isotropic Friedmann- Robertson- Walker 
universe given by the line element 

ds 2 = -dt 2 + a 2 (t)S l:j dx l dx j , (3) 
where a(t) is the scale factor. The Friedmann equation and the continuity equation read 

3H 2 = E, (4) 

E=-3H(E + P), (5) 

where dot denotes derivative with respect to cosmic time, the Hubble rate is H = a /a, E is the energy of the inflaton 
and it is given by 

E = 2XP X - P, (6) 
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where P x denotes the derivative of P with respect to X. It was shown in [7(| (see also [73) that for this model the 
speed of propagation of scalar perturbations ( "speed of sound" ) is c s given by 

c * = ^ = (7 ) 

s E yX P x + 2XP XX • V ' 

We define the slow variation parameters, analogues of the slow-roll parameters, as: 

H XP X e c s 

e = -^ = ^' ? = ^r S = z^- (8) 

We should note that these slow variation parameters are more general than the usual slow-roll parameters and that 
the smallness of these parameters does not imply that the field in rolling slowly. We assume that the rate of change 
of the speed of sound is small (as described by s) but c s is otherwise free to change between zero and one. 
The power spectrum of the primordial quantum fluctuation was first derived in [70| and reads 

1 E 2 1 H 2 

where it should be evaluated at the time of horizon crossing c ssr k = a^H*. The spectral index is 



_ dhxP c (k) 



dink - 2£ -^ S - < 10 > 

WMAP observations of the perturbations in the CMB tell us that the previous power spectrum is almost scale 
invariant therefore implying that the three slow variation parameters should be small at horizon crossing, roughly of 
order 10~ 2 . 



III. NON-LINEAR PERTURBATIONS 



The main goal of this section is to calculate the total connected trispectrum of the primordial curvature perturbation 
for the model ([1]), at leading order in the slow-roll expansion and in the limit of small sound speed. In the next 
subsection, with very simple arguments, we will estimate the order of magnitude of the trispectrum coming from the 
scalar exchange and compare it with the result of the trispectrum coming from the contact interaction [25| | (see also 
[30|). In subsection IIII B[ we review the calculation of the contact interaction trispectrum. Finally, in subsection 
IIII C[ we calculate the scalar exchange trispectrum. The total trispectrum will be the sum of these contributions. 

Non- linear perturbations of the model |T]) have been considered by several authors before. In [14| . Gruzinov obtained 
the bispectrum using the simple method of expanding the action of the inflaton only. This method gives the leading 
order bispectrum in slow-roll and in the small sound speed limit. In [24[, the authors performed a detailed analysis 
of the bispectrum produced in the present model including also next order terms in slow-roll and small sound s pee d 
expansion. The trispectrum has also been studied in the literature. In particular, Chen et al. [25j (see also [301 ]) 
calculated the trispectrum from the contact interaction. 

The Lagrangian density of the inflaton expanded up to fourth order in perturbations (at leading order in slow-roll 
and in the small sound speed limit) can be founded in [25| . and it reads 

£ (2) = ^ 2 -^(**) 2 , (11) 



A3 



£ (3) = Pxx^ + Pxxxy « 3 " 3 - P,xx^aa(8a) 2 , (12) 



,(4) _ £ 
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?-P AX + 4> Z P, XXX + -P XX a 4 - a- A 6t\day (j) 2 P xxx + P XX + -a^P xx (da 



6 " T ~ '"' w * 2 

where a denotes the field perturbation 8<fi(t, x) 
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(13) 
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In the next subsections, we will use the so-called "interaction picture formalism" [73[ to calculate the four-point 
quantum correlation function. In this formalism one needs to know the interaction Hamiltonian. The Hamiltonian 
density H can be found using Eqs. (fTTjl. (fT2| . (fT3| and 

H = nd - C, (14) 

where C is the Lagrangian density and n is the momentum density, defined by 7r = dC/dd. The third order Hamil- 
tonian density is Ti^ = — £' 3 ' and the fourth order Hamiltonian density is [25| 

U {i) = /3ia 3 d 4 + f3 2 aa 2 (da) 2 + faa' 1 {daf , (15) 

where we have used d = dTlo/dir to express 7r in terms of d, and Tig denotes the kinematic Hamiltonian density that 
is quadratic in tt and a. Finally to obtain the interaction Hamiltonian from the (|15j) . the variables d and a should 
be replaced with their interaction picture counterparts di and a/. 
The constants are defined as 

ft = ~P,xx (l - \c 2 )j + \j> 2 c 2 s P X xx, 

ft - -| Pxx- (16) 
For DBI inflation, at leading order in the sound speed, the constants are 

2c^0 2 4cf</> 2 8c s 2 

From the second-order action, we can solve for the perturbation and quantize it according to the standard procedures 
of quantum field theory: 

afa.x) = ^3 J d 3 k[u( V ,k)a(k)+u*( V ,-k)a^~k)}e^ , (18) 

where u(r),k) = —j=—p= — (1 + i/cc s ?7)e~ 4fcCs '' is the solution of the quadratic Lagrangian, and r\ = —-^jj is the 
conformal time. The commutation relations are given by 

[o(ki),a(k 2 )] - [a t (k 1 ),a+(k 2 )] = 0, [a(k x ), a t(k 2 )] = (2vr) 3 ^ 3 )(k 1 - k 2 ). (19) 



A. Estimate of the amplitude of the total trispectrum 



In this subsection, we shall present a very simple argument to prove that the trispectrum coming from the contact 
interaction diagram is of the same order as the trispectrum coming from the scalar exchange diagram. Therefore the 
total trispectrum is the sum of these two contributions. For simplicity, we will consider the DBI-inflation case but 
the conclusion is valid for generic models of the form (fTJ). 

Let us start by drawing the tree-level diagram that contributes to the bispectrum. It can be found in Fig. [TJ and 
it is composed of three external lines labeled by C and it has three scalar propagators Vq ~ H 2 e _1 cj 1 , and one third 
order vertex I3 as depicted in the l.h.s. of Fig. [5] The bispectrum coming from this diagram is then ~ V^h- 
The third order vertex is given by ^3 ~ H~ 2 tc~ x as is easily seen from the third order Hamiltonian [24|, thus the 
DBI-bispectrum scales like Bq ~ c~ 2r P 2 (or equivalently the non-linearity parameter Jml scales like /jvl ~ c ^ 2 )- F° r 
the trispectrum, there are two relevant tree-level diagrams, Fig. [3J If we denote the fourth order vertex on the r.h.s. 
of Fig. H]by I4 then the diagram on the l.h.s. of Fig. Ogives a contribution for the trispectrum as T CI ~ V^I^ and 
the diagram on the r.h.s. of the same figure gives T SE ~ V^I^, From the fourth order Hamiltonian, the fourth order 
vertex is estimated as I4 ~ H~ 2 ecJ 3 thus the trispectrum from the contact interaction is given by T CI ~ cJ 4 V 3 [25| . 
The magnitude of T SE is easily calculated as T SE ~ cJ 4 V 3 which has the same magnitude as T CI . Therefore for 
consistency the total trispectrum is the sum of T GI and T SE . 
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FIG. 1: Bispectrum diagram 




Recently, [6!| showed that in the case of the standard kinetic term inflation, i.e. when P is P = X — V((f>) (V(4>) 
is the potential), the diagram where we have an exchange of a graviton can in fact dominate the contact interaction 
diagram. In this model, the contribution from the scalar exchange diagram is negligible in comparison with both the 
contact interaction and the graviton exchange. They showed that the contact interaction trispectrum is T CI ~ eV^, 
the graviton exchange trispectrum is T GE ~ eP^ and the scalar exchange trispectrum is T SE ~ e 2 V^. The third 
order scalar vertex is 7^3 ~ H~ 2 e 2 and the fourth order vertex is 1^4 ~ H~ 2 e 2 . 

We have checked that in the case of non-standard kinetic terms the graviton exchange diagram gives a contribution 
to the trispectrum that is suppressed by e and c s with respect to the leading order and we neglect this contribution. 
This fact can be easily understood if one uses the main result of [69( and one recalls Maldacena's result [4[ that the 
third order vertex of two-scalars-one-graviton is I 1 ^2 ~ H~ 2 ec s . 

In the next subsections, we will derive the momentum dependence of the trispectrum from the contact interaction 
and the scalar exchange. 



B. The trispectrum from the contact interaction 

In this subsection, we shall make use of the so-called "interaction picture formalism" [73[ to calculate the four-point 
function of the field perturbation at horizon crossing coming from a diagram like the one on the l.h.s. of Fig. [3] The 
result can also be found in Chen et al. [25[ (see also |3(J). These results are at leading order in slow-roll and in the 
limit of small sound speed. For the contact interaction diagram, the four-point function is 



( W(0, ki)ty(0, k a )ty(0, k 3 W(0, k*)|n) 



ci 



dr,(Q\ [<5</» / (0,k 1 )^ 7 (0,k 2 )^ / (0,k 3 )^ / (0,k4),i7j 4) (?7)l |0), (20) 



where Hi{rj)^' is the fourth order interaction Hamiltonian given can be found from Eq. |15 

and reads 



(21) 



where prime denotes derivative with respect to conformal time rj. 



FIG. 3: Trispectrum diagrams. On the left: the contact interaction. On the right: The interaction via exchange of a scalar. 



For a generic K-inflation model, we obtain 

(n|ty(o, ki)ty(o, k 2 )^(o, ka)ty(o, k 4 )\n) CI 

1152/3i7V 8 c^ 



-(27r) 3 <5^(K) 



1 



fc2fc 2 2 (k 3 -k 4 ) 



E 



12- 



fc 3 fc 4 



+32/3; 



TV 8 (k 1 -k 2 )(k 3 -k 4 ) + (k 1 .k 3 )(k 2 -k4 



(EL** 
(ki -k 4 )(k 2 



permutations 



k 3 ) 



Ei=i fe» nf =1 fc 3 



V 



(Etr *i) (Ett * 



2 — 1 



- + 12- 



(22) 



where the total momentum K is defined by K = k 4 + k 2 + k 3 + k 4 , the constant N is defined by N — Hj yj2Pxc s 
and "permutations" denote the other twenty three terms that result from the permutations of {fci, fc 2 , fc 3 , /c 4 } in the 
preceding term. The three dimensional Dirac delta function of K ensures momentum conservation and implies that 
the four momentum vectors k^, where i = 1, . . . , 4, form a quadrilateral. 

For DBI-inflation, in the small sound speed limit, the terms proportional to /3 2 and /? 3 are sub-leading with respect 
to the term proportional to f3\. For K-inflation models, this is not necessarily the case and this fact can be used to 
distinguish some K-inflation models from DBI-inflation [25j ] . 

To obtain the four-point correlation function of the curvature perturbation at some time after horizon crossing we 
can use the linear relation Q = —(H/(f>)S(f> because the higher order terms in this relation only generate sub- leading 
corrections to this result and we ignore them. Finally, we obtain 



H 4 



(0|C(o, kOC(o, k 2 )c(o, k 3 )c(o, k 4 )|^) 0i = -^( W(o, ki)<ty(o, k 2 )^(o, k 3 )^(o, k 4 )|Q> 



CI 



(23) 



C. The trispectrum from the scalar exchange interaction 

Within the "interaction picture" formalism, to calculate the four-point function resulting from a correlation estab- 
lished via the exchange of a scalar mode as depicted in the r.h.s. of Fig. [3] one needs to evaluate the following time 
integrals 



(fi|ty(0, ki)(ty(0, k 2 )<J0(O, k 3 )50(O, k 4 )|f>) 



SE 



— oo •/ — oo 



dr, / dfj{0\ \ \SM0, k 1 )^(0,k 2 )^ J (0,k 3 )^(0,k 4 ), J fff ) (77)1 ,H^(fjj\ |0), 



(24) 



where Hj 3 \r]) is the third order interaction Hamiltonian given by Hj 3 ^ — J d^xH^p — — J d 3 xC^p . We will use 



Hf\r 1 )= I d 3 x AaScpf + Ba6$ (dSfaf 



(25) 
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where the constants A and B are defined as 

1 



A = -"-lPxx + ^ 2 Pxxx), B = ^P XX . (26) 



For DBI inflation they are given by 



A= !— , B=-^—. (27) 

20c* 20cf 



The calculation of the integral in (JM|) is rather long but relatively straightforward (however see uM)- The total 
four-point function of the field perturbations is 

2 AT 4 

(Q|^(o,k 1 )^(o,k 2 )^(o,k 3 )^(o,k 4 )|n) S£; = {2k)H^(y: ~ 



(kik 2 k 3 k 4 )2 

9A 2 ( ^i(fci,fc2, —ku, k 3 , k 4 , k 42 ) - Ti(-ki,-k 2 , -km, k 3 , fe 4 , ku) 



+AB ^3(k 3 ■ k 4 ) \^F 3 (h,k 2 , -h 2 ,ki 2 ,k 3 ,k 4 ) - F 3 (~k 1 ,-k 2 , -k 12 ,k 12l k 3 , fc 4 

+6(ki2 • k 4 ) ^V 3 (fc l7 fc 2 , -fci2, fc 3 , fe 4 , feia) - T 3 (-ki, -k 2 , -k 12 , k 3 , fc 4 , fei 2 
+3(ki ■ k 2 ) f.F 4 (-fci 2 , h,k 2 ,k 3 , ki, k 12 ) - T 4 (-k 12 ,-k 4 , -k 2 , k 3 , k 4l k 12 ) 
-6(ki2 ■ k 2 ) KF 4 (fei, k 2 , -k 12 , fc 3 , k 4 , k 12 ) - F 4 (-ki, -k 2 , -k 12 , fc 3 , fc 4 , fei 2 

-B 2 ^(ki • k 2 )(k 3 • k 4 ) ^F 2 (-ki 2 , ki, k 2 , k\ 2 , k 3 , k 4 ) - F 2 {-ki 2 , —ki, -k 2 , k X2 , k 3 , fe 4 ) 

+2(ki • k 2 )(ki2 ■ k 4 ) \^F 2 (-ki 2 , fei, k 2 , k 3 , fe 4 , fc i2 ) - F 2 (-k 12 , -k x , -k 2 , k 3 , k 4 , k 12 ) 
-2(ki2 • k 2 )(k 3 • k 4 ) I T 2 (ki,k 2 , -k 12 , k 12 , k 3 , k 4 ) - F 2 (-ki,-k 2 , -k 12 , k 12 , k 3 , k 4 ) 



-4(ki2 • k 2 )(ki2 • k 4 )\^F 2 (ki,k 2 , -k 12l k 3l k 4l k 12 ) - F 2 (-ki, -k 2 , -k 12l k 3l fe 4 , fc i2 ) 
+23 permutations of {k\, k 2 , k 3 , k 4 }, (28) 



where as before K = ki + k2 + k 3 + k 4 , N — Hj ^/2Pxc s and k\ 2 = |ki2 1 = | ki + k2 1 . Similarly, we define k a i, as 
k a b = |k Q b| = |k a + k(,|, where k a and kh represent any of the four momentum vectors ki, k2, k 3 and k 4 . Momentum 
conservation implies fei2 = k 34 , k\ 3 = k 24 and fci 4 = k 23 . 

Because there are only four different types of double time integrals over the mode functions, we define four Ti 
functions (with i = 1, . . . , 4). Their explicit expressions are 

Fi{k u k 2 ,k 3 M,kM = / d-qa^) /' d^a{ii)U*\i l M)U*\vMU*\riM)U*\iiM)U'\nM)U'\nM 



*6|*i fl + 3^+6^]. (20, 



H 2 U| A 3 C 3 V C C 2 

where A is defined by the sum of the last three arguments of the Ti functions as A = k 4 + k 5 + k 6 and C is defined 
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by the sum of all the arguments as C = k% + k 2 + k 3 + k 4 + &5 + kg- The remaining functions are 



Ajofo) / df,a{fi)U* (r), ki)U*(<n, k 2 )U*{ V , k 3 )U* (fj, h)U*(fj, fo)U*(fj, fc 6 ) 



oo J — oo 



ki fc 4 | 2 1 



H 2 \k 2 k 3 k 5 k 6 \l AC 



1 + + 2- 



A 



A 2 



+\ (h + k 3 + k 5 + k 6 + 1 ((fe 2 + fe 3 ) (*8 + fee) + 2fc 5 fc 6 ) + 2 k5k6{k ^ 2 +h) 



+g2 yk 5 k 6 + (k 2 + k 3 ) (k 5 + fc 6 ) + k 2 k 3 

+ \ (k 2 k 3 (k 5 + k 6 ) + 2k 5 k e (k 2 + k 3 )) + 2 k2hk f 6 



+ ^ (k 2 k 3 (k 5 + k (l ) + k 5 k 6 (k 2 + k 3 ) + 2 k2hk5k6 



24 



k 2 k 3 k 5 k 6 
C 3 



(30) 



T 3 {ki,k 2 ,k 3 ,k 4 ,k 5 ,ko) 



J r 4{ki,k 2 ,k 3 ,k 4 ,,k 5 ,k 6 ) 



dr,a( v ) f V dfja{fj)U*' { v , kx)U*' { v , k 2 )U*' (r?, k 3 )U*' [fj, k 4 )U*(fj, h)U*{fj, k 6 ) 

-oo J —oo 

2 — — + — I k 5 + k e + 2 



H 2 \k 5 k 6 \i AC S 



1 



A 



A 2 C 



A 



12- 



k 5 k 6 



C 2 



(31) 



d V a( V ) P dfja{fj)U*' (v, fci)t/*(?7, k 2 )U*( V , k 3 )U*' {fj, k 4 )U*' (fj, k 5 )U*' (fj, k e ) 



= 2 



N 6 c 4 s \hkihke\i 1 



H 2 



\k 2 k 3 



A 3 C 



A A?_ k 2 + k 3 n A(k 2 + k 3 ) + k 2 k 3 
+ c + c2 + c + c2 

+3^ (A (k 2 + k 3 ) + 2k 2 k 3 ) + 12fc 2 fc 3 ^ 



where the function U(rj,k) is a modified mode function given by 



(32) 



U(r), k) = N— (1 + ikc s T]) e~ lkc ^. 
IfcU 



(33) 



If the sign of the argument k of U is positive then U is equal to the mode function, if the sign is negative then U 
equals the complex conjugate of the mode function. 

The study of the momentum dependence of the trispectrum will be the subject of the next section, but here we would 
like to make a remark. In the configuration where all the momentum vectors have equal magnitude, it is easy to see 
that the variable C defined previously and that appears in some Ti functions in (f28l) becomes zero. Because C appears 
in the denominator of some fractions in the definitions of Ti one might naively think that the trispectrum diverges 
for these configurations. This is not the case. In Appendix 13 we show that if one includes all the permutations 
of these terms the apparent divergences do not appear anymore and as expected the trispectrum is finite for these 
configurations. 

As in the previous subsection, the four-point function of the curvature perturbation £ is related with the four-point 
function of the field perturbation as 



<fi|C(o, ki)c(o, k 2 )c(o, k 3 )C(o, k 4 )|ft) SB = ^-(n\s<t>(o, ki)^(o, k 2 )^(o, k 3 )^(o, k 4 )|ft> 



SE 



(34) 



This constitutes one of the main results of this work. 

The total trispectrum of Q is the sum of the contributions coming from both of the diagrams of Fig. [3] as 

(o|C(o, ki)c(o, k 2 )C(o, k 3 )c(o, k 4 )\n) total = (n|c(o, ki)c(o, k 2 )C(o, k 3 )c(o, k 4 )\n) CI 

+ (O|C(0, ki)C(0, k 2 )C(0, ka)C(0, k 4 )|fi) SE . 



(35) 
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From the previous result one can read the order of magnitude of the total trispectrum and it agrees with the 
discussion of subsection IIII Al Moreover, one can also read the exact momentum dependence of the trispectrum. 
In the next section we shall study the shape of the trispectrum in more detail. In particular, we will consider the 
so-called equilateral configuration and we shall calculate the value of the non-linearity parameter tnl- 



IV. tnl AND THE SHAPE OF THE TRISPECTRUM IN THE EQUILATERAL CONFIGURATION 

In this section, we will calculate the non-linearity parameter tjvl i n the limit where all the four momentum vectors 
have the same magnitude k. This is the so-called equilateral configuration used in previous works [25l . l68j . If we 
denote the angle between k^ and kj by 0y, then in this configuration, we have cos(#i2) = cos(#34) = cos(0 3 ), 
cos(023) = cos(6*i4) = cos(6*i), cos(#i3) = cos(6*24) = cos(#2) and cos(#i) + cos(#2) + cos(^a) = —1 due to momentum 
conservation. 

The inflationary trispectrum is often parameterized [l^, H^, [6^| by two non-linearity parameters T l § c L al and g l § c L al as 



(C(k 1 )C(k 2 )C(k 3 )C(k 4 )) = T%l al (T c (k 13 )T c (k 3 )V ( (h) + 11 permutations) 

+ %9 l NL l (M^W^Vcih) + 3 permutations) , (36) 

where the power spectrum VAk) is given by (C(ki)C(ka)) = (2Tr) 3 S^(k 1 +k 2 )V c (k 1 ) or V c {k) = H 2 /(4ec s k 3 ). This 
shape of the trispectrum results from the most general local parametrization for the non-linearity of £ as 

C(x) = C G (x) + \{t^Y /2 (CS(X) - CSM) + ^9w?& «> ( 37 ) 



where Cg( x ) denotes the first order curvature perturbation (gaussian random variable), Cg( x ) denotes the average of 
(^q. The trispectrum calculated in the present work results from quantum correlations around horizon crossing and 
does not admit a simple parametrization in terms of momentum independent non-linearity parameters as (|36p . 
In the next subsections, we will adopt a naive parametrization for the trispectrum as 

r J vL(ki,k2,k3,k 4 ) = f^r) T c (ki,k 2 ,k 3 ,k 4 )Ht = i£:.? 



{k\k2 + ^3^4) {k\3 + ^14 3 ) (^1^4 + ^2^f) (^12 3 + ^13 3 ) + (^1^3 + ^2^1) (^12 3 + ^14 3 ) 



-1 



(38) 



where Tf is related to the trispectrum of C, as 

(C(k 1 )C(k 2 )C(k 3 )C(k 4 )) = (2 7 r) 3 ^ 3 )(K)T c (k 1 ,k2,k 3 ,k 4 ). (39) 

In general, tml defined by the previous expression will be shape dependent. For instance, this will be the case for the 
contact interaction trispectrum. Only for non-Gaussianity of the form ([57)) with g l ^ L al — the parameter tjvl will be 
a constant, i.e. shape independent. If one fully specifies the momentum vectors configuration, then tnl becomes just 
a number. This seems to be the approach followed by several authors in the literature and in the next subsections we 
shall do the same to obtain the order of magnitude of the trispectrum. However, we should emphasize that the full 
trispectrum contains much more information than the one that can be described by just a number like tjvl and that 
this parametrization might not be so adequate to describe the information encoded in a general trispectrum. 



A. tnl- general model 

For the general model (pj and for an angular configuration like cos(#i) = cos(#2) = cos(#3) = —1/3 the contact 
interaction trispectrum gives a non-linearity parameter as 

t N l% ~ -0.14^J^/3i, t nl % ~ 0.2l|//? 2 , t nl % ~ -0.55^^/33, (40) 

where the different values of tnl correspond to the three terms in the contact interaction trispectrum Eq. 
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The scalar exchange trispectrum gives values for t^l as 

tnl S A 2 ~ 0.52^^ A 2 , t nl s a e b ~ -1.75^J^AB, r WL |f ~ 2.51^B 2 , (41) 

,VC 

where rjvi^f , tnl b % and Tjvx^ correspond to the values of tjvx calculated from the terms in (|28p proportional to 
A 2 , B 2 and AB respectively. 



B. The shape of the trispectrum for the DBI-inflation model in the equilateral configuration 

In this subsection, we shall study the shape of the trispectrum for the DBI-inflation model in the equilateral 
configuration. We will assume that the trispectrum is maximized in the equilateral configuration and we will study 
the dependence of the maximum on the two remaining degrees of freedom, the angles Q\ and O2 1 ■ We will show that 
the trispectrum is maximized in the configuration where the angles between the four different momentum vector are 
equal. Although we do not prove that the trispectrum has the maximum value for the equilateral configuration, we 
present some evidence indicating that this might be the case in Appendix [Bl 

In the equilateral configuration, the contact interaction trispectrum is constant and only the scalar exchange 
trispectrum depends on 61 and (9 2 . The plot of the shape of the scalar exchange trispectrum can be found in Fig. 
HI The exact analytical expression can be found in Appendix [Cj The shape of the trispectrum has a maximum for 
cos#i = cos 6 2 — cos 6*3 = —1/3 as it can be better seen in the density plot of the trispectrum in the r.h.s panel of 
Fig. [D In Fig. we plot the non-linearity parameter rf^ calculated from the scalar exchange trispectrum. The 
maximum amplitude occurs for cos#i = cos #2 = cos #3 = —1/3. 

When cos 9% = — 1 or cos 02 = — 1 or cos #3 = — 1 (on the diagonal line in the plot of Fig. [5l when cos #2 = — cos#i) 
the non-linearity parameter is zero. This result can be nicely understood using the counter-collinear limit of the scalar 
exchange trispectrum (o9j . According to [6^ | , in the limit where the momentum of the scalar mode that is exchanged 
goes to zero, one can find a simple relation between the scalar exchange trispectrum and the power spectrum, in 
a similar way to Maldacena's consistency relations [1] (see also [1^, ^M, [ZSj]) since one can treat this mode as a 
background. Summarizing their result, let us suppose that the momentum of the exchanged particle is k\2 and that 
ki2 <C fci « fc2,^3 ~ fed- Then the mode associated with this scalar particle will cross the horizon much before the 
other hi modes, where i = 1, . . . , 4, and its only effect is to rescale the spatial background where the ki modes exist. 
After some algebra, Ref. [69( found that the scalar exchange trispectrum obeys the following relation 

<C(k 1 )C(k 2 )C(k 3 )C(k4)) SB - (2^)V 3 )(K)(n s - l) 2 7\(fci2)n(W(fc3), (42) 

in the limit ki2 — > 0. The previous equation implies that in the counter-collinear limit, the scalar exchange non- 
linearity parameter Tjyf is of order e 2 (at leading order in our approximations this is equivalent to say that rjyf 
should vanish). We have verified that Tjyf calculated from (|2"5)l does indeed vanish in the limit k a i> — > 0. This 
provides a consistency check on (|28l) . If we note that in the equilateral configuration ^14 = ^23 = \/2k^/l + cos6*i, 
k 2 4 = fci3 = \/2fc\/l + cos 6*2 and /C34 = k i2 = \[2k\/\ + cos6* 3 , then it is easy to understand why in the plot of Fig. 
t nl vanishes for cos#j = —1 with j = 1, 2, 3. This is because we are in the domain of validity of the counter-collinear 
limit relation. 



C. tnl: DBI-inflation 



In order to give numerical results for tjvl j one has to choose a particular model which will determine the parameters 
(3i, A and B. In this subsection, we particularize the results (|4"0")) and (|4T|) for the DBI-inflation model at leading 
order in the small sound speed. 

The total tnl coming from the scalar exchange trispectrum is 

SE ^ SE,^ SE,^ SE 0- D fA o\ 
T NLDBI — T NL A 2 + T NLB 2 + TNLAB ~ — j - ' ^1 



1 Contrary to what happens in the case of the bispectrum, the equilateral configuration conditions do not fix all degrees of freedom (dof) 
required to describe the shape of the trispectrum and we are left with two angular dof [88] 
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FIG. 4: Left: The shape of the scalar exchange trispectrum as a function of the variables cos#i and cos 62- The maximum 
amplitude has been normalized to unity. Right: Density plot of the shape of the scalar exchange trispectrum as a function of 
the variables cos#i and cos #2- The maximum amplitude occurs for cos Si = cos 62 = —1/3. 




FIG. 5: Plot of the non-linearity parameter calculated from the scalar exchange trispectrum. The maximum amplitude 
occurs for cos Si = cos #2 = — 1/3 and its value was rescaled by c 4 a . 



This result should be summed with tjvl coming from the contact interaction trispectrum [25| 

-r ci ci 0.04 

TNLDBI = T NLf3 ± ~ —■ (44) 

c s 

For this particular configuration the magnitude of tnl^xbi is more than one order of magnitude higher than tnlz/bI' 
We conclude that the total non-linearity parameter rjv l for DBI-inflation in the equilateral configuration with equal 
angles between the momentum vectors is 

n ^fi 

tnl£ bi - —r- (45) 



V. CONCLUSION 



In this work we have computed the connected four-point function of the primordial curvature perturbation at 
leading order in the slow-roll expansion and in the small sound speed limit for a generic model of kinetically driven 
inflation. This four-point correlation is coming from either via the exchange of a scalar mode or by a four-leg contact 
interaction at a point. We showed that in general the scalar exchange correlation is as important as the contact 
interaction correlation previously calculated in the literature. So the final answer to be compared with observational 
data should include the sum of these two contributions. 
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We have studied the shape dependence of the trispectrum for the DBI-inflation model in the equilateral configu- 
rations. The equilateral limit conditions do not fix completely the configuration of the four momentum vectors and 
we are left with two remaining angular degrees of freedom. We have shown that the trispectrum is maximized when 
the angles between the four different momentum vectors are equal to arccos(— 1/3) w 109.5°. We found that for 
this "maximal" equilateral configuration the non-linearity parameter for DBI-inflation tni^dbi * s arj out one order 
of magnitude larger than the existing result in the literature that only takes into account the contact interaction 
diagram. Our final result for the non-linearity parameter for the DBI-inflation model is TjVi^xBj ~ 0.56/c^. 

We did not prove that the maximum of the trispectrum occurs for the equilateral configuration. But we gave 
some evidence in Appendix B. In general there are five parameters to characterize the shape of the trispectrum. It 
would be desirable to explore, possibly numerically, all configurations to show that the maximum of t^l lies in the 
configuration we found. Also it is necessary to construct an estimator to extract information from CMB anisotropies. 
We will leave these issues for future investigations. 

Note added: On the day this work was submitted, the paper (77j appeared in the arXiv, which contains some similar 
results. 



Acknowledgments 

We would like to thank Filippo Vernizzi for interesting discussions. We thank Takashi Hiramatsu for pointing 
out some typos in an early version of the manuscript. FA and SM are supported by the Japanese Society for the 
Promotion of Science (JSPS). KK is supported by RCUK, STFC and ERC. The authors thank the Yukawa Institute 
for Theoretical Physics (YITP) at Kyoto University. Discussions during the YITP workshop YITP-W-09-01 on 
"Non- linear cosmological perturbations" were useful to complete this work. KK's visit to the Yukawa Institute was 
supported by the Royal Society international joint project grant. We also thank the Institute for the Physics and 
Mathematics of the Universe (IPMU) for organizing the workshop entitled "Focus week on non-Gaussianities in the 
sky" during which we had several important discussions about this work. SM is grateful to the ICG, Portsmouth for 
their hospitality when this work was completed. 



APPENDIX A: CANCELING THE APPARENT DIVERGENCES OF THE TRISPECTRUM IN THE 

EQUILATERAL LIMIT 

In the configuration where all the momentum vectors have equal magnitude, i.e. the equilateral configuration, the 
variable C defined as the sum of the arguments of the T functions vanishes for the terms in (|28p that contain a Ti 
function with three negative arguments. Because C appears in the denominator of some fractions in the definitions of 
T% one might naively think that the trispectrum (|28|) diverges for these configurations. We will now show that when 
considering all the 23 permutations in ([25]). the powers of C that appear in the denominator in the different terms 
cancel out and we are left with expressions that are clearly finite in the limit C going to zero. 

The proof is rather simple if one notes that the apparently divergent terms (Ti functions with three negative 
arguments) can be gathered in certain combinations that are clearly finite in the equilateral limit. 

For example, if we take all the divergent terms containing Ti, we find pairs of terms like 
T\(— k a , — kb, — k a b, k c , kd, k a b) + Ti(—k c , — kd, ~k a b, k a , kb, k a b), where the only difference from the first term to the 
second is that k a ,kb changed position with k c ,kd, where k a ,kb,k c ,kd can represent any of the momentum vectors 
ki, k%, k$, &4 and k a b is k a b — |k a &| = |k a + k&|. To complete the proof, one just needs to show that this combination 
is finite in the equilateral limit. After some algebra one can show that 

N 6 c® 2 l 1 

T\(— ka-, —kb, —k a b, k C: kd, k a b) + T\(—k c , kd, ~k a b, k a , kb, k a b) — ^—^2-\k a kbk c kdk ab \ 5 ^3^3 , (Al) 

where A\ = —(k a b + k c + kd) and A2 = — (k a b + k a + kb). The r.h.s. of the previous equation is clearly finite in the 
limit C -> 0. 

The proof of finiteness for the remaining divergent terms containing T3 and Ti in (|28p is completely analogous to 
the proof just described but now one has to use the following combinations of terms. 

Ts( k a , kb, k a b, k a b, k c , kd) + T^( k a b, k c , kd, k a , kb, k a b) 

= 2 N e cj \k a k b kl b \h 

H 2 \k c k d \ % 



k c kd fee "t~ kd 



J A 3 A 3 



A 2 A 3 



(A2) 
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•F~s( ka, k b , k ab , k c , kd, k ab ) + 3~a\ k c , kd, k ab , k a , kb, k a b) 



,A 6 cf \k a k b k c k ab \2 
' H 2 \k d k ab \i 



1 



kdk. 



d K ab 



A, A3 A3 A3 1 A2 A3 



(A3) 



Finally, all divergent terms containing Ti can be completely written in terms of the following combinations 



•F~i( k ab , k a , k b , k ab , k c , kd) + k ab , k c , kd, k ab , k a , k b ) 

\k ab \ [2A\ + A 2 (k ah - 2k a ) - 2k a (k a + kab)] [2A\ + Ax (k ab - 2k d ) - 2k d (k d + k ab ) 



iV 6 c? 



H 2 \k a k b k c k d \^ 



A 3 A 3 



(A4) 



•F~i( k ab , k a , k b , k ci kd, k ao ) + J~%( k c , kd, k ab ,k ab ,k a ,k b ) 



\k c k a b 



[2AI + A-2 (k ab - 2k a ) - 2k a (k a + k a b)] [A\ - Al (k d + k ab ) + 2k d k ab ] 



H 2 \k a kbkdkab\i 



A 3 A 3 



(A5) 



J~2( k b , k a , k ab , kd, k c , k ab ) + J~%( kd, k c , k ab , k b , k a , k a b) 

N e c 2 [fcbfcrflj [2A\ + Ax (kg ~ 2k ab ) - 2k ab (k d + k ab )] [A\ - A 2 (k a + k ab ) + 2kgk ab ] 
H 2 \k a k c k 2 ab \i A\Al ■ {Ab) 

All of these combinations are finite in the equilateral limit and all apparently divergent terms in the scalar exchange 
trispectrum (J2HJ) can be written in terms of these combinations. This implies that the trispectrum is also finite 
in this limit as expected. 



APPENDIX B: WHERE DOES THE MAXIMUM OF THE TRISPECTRUM FOR THE DBI-INFLATION 

MODEL LIE? 

It seems difficult to pursue an analytical answer for this question. The main reason is the fact that we need five 
degrees of freedom to parameterize the shape of the trispectrum and the expression for the trispectrum, Eq. (|55|l . is 
fairly complicated. Numerical searches for the maximum in the full parameter space are possible even with such a 
large parameter space to explore. We leave this for future work. 

In this Appendix, we will relax the condition of equilateral configuration, allowing for two of the momentum vectors 
to have a different magnitude from the other two and we will show that the maximum that we found in subsection 
IIVBI is still a maximum of tjvl • This provides some evidence that this particular configuration maximizes l ■ 

Before we show the results, it is interesting to consider analytically certain limiting configurations such as the 
squeezed limit or in the case of the scalar exchange diagram the counter-collinear limit. These limits give us other 
indications about where the maximum of the trispectrum might lie. Also, because these limits can be obtained using 
different methods from the one used in this work, they provide consistency checks with our final answers. 

As first pointed out by Maldacena [3| and Seery et al. [68| . all higher order correlators of £ in single field inflation 
obey consistency relations in the squeezed limit, i.e. when one of the momentum vector is very small. Let us consider 
the limit ki — > 0, which implies that the corresponding mode C(ki) leaves the horizon well before all other modes 
leave the horizon. By the time the remaining modes exit the horizon, C(ki) will be frozen as a super- horizon mode 
and its only effect is to deform the background. Following the algebra with some more detail, one finds the following 
relation for the trispectrum [68| 

(C(k 1 )C(k 2 )C(k 3 )C(k 4 )> - -p c (fc 1 )-^(C(k 2 )C(k 3 )C(k 4 )>, (Bi) 

in the limit ki — > 0. Using the fact that the bispectrum is Bq ~ cJ 2r P 2 it follows that in the squeezed limit according 
to (IB1|) the trispectrum is of 0(ec~ 2 V^). Comparing this result with the full trispectrum, Eq. which is of order 

e -3 , one finds that the total non-linearity parameter t^l calculated from (|35|) should be at most of order ecj 2 (i.e. 
under our approximation it should vanish) when any of the momentum vectors goes to zero. We verified that this is 
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the case for (I35I) 2 . This would mean that the maximum of tjvl might lie near the equilateral configuration. 

As an example, we considered configurations where the two momentum vectors IC3 and IC4 have a different magnitude 
q from the magnitude of ki and k.2 that we normalized to unity. Then in this configuration, we have cos(#34) = 
cos($3), cos(6>23) = cos(#i4) = cos(6*i), cos(#i3) = cos(#24) = cos(6 2 ), cos(#i2) = q 2 (1 + cos(03)) — 1. The momentum 
conservation implies cos#3 = — (l/<7)(cos#i + cos#2) — 1- We also allow one of the angular dof to vary freely, the other 
angular dof 9\ is fixed by cos#i = —1/3 to the value that gives us the maximum that we found in subsection IIV Bl 
In Fig. [5J we plot the non-linearity parameter tjvl for the DBI-inflation model. From the momentum conservation, 
cos(#2) is restricted to cos(#2) < 1/3. Also we restricted q as q > 1 as the configurations with q < 1 can be related 
to those with q > 1 by rescaling k\ and k%. The left panel shows the contribution from the contact interaction. In 
the right panel, we plot the total tjvl- As can be seen in the r.h.s panel of Fig. the maximum of Tj^°^ al occurs 
exactly for the equilateral configuration when all the angles 9\ , 9 2 and 9% are equal. This provides evidence that the 
maximum of the non-linearity parameter would happen for this configuration. 



APPENDIX C: THE DBI SCALAR EXCHANGE TRISPECTRUM IN THE EQUILATERAL LIMIT 

In this Appendix, we present the scalar exchange trispectrum of the primordial curvature perturbation Q for the 
DBI-inflation model in the equilateral configuration. These analytical expressions can be derived from the general 
expression for the scalar exchange trispectrum, Eq. (f2"5)) . We will use these simpler equations to calculate the non- 
linearity parameter tnl in subsection IIV CI and to study the shape dependence of the trispectrum in subsection 

WM 

At leading order in the small sound speed and in the slow-roll expansion, the scalar exchange trispectrum for the 
DBI inflationary model is 

(fi|C(o, k x )c(o, k 2 )C(o, k 3 )C(o, k 4 )\n) SE = (O|c(o, k x )c(o, k 2 )C(o, k 3 )C(o, \^m SE A2 

+<ft|C(o, k x )c(o, k 2 )C(o, k 3 )C(o, k 4 )\n) SE B2 
+{n\ao, kOCCo, k 2 )C(o, k 3 )C(o, k 4 )\n) SE AB , (ci) 

where the three different contributions come from the terms in that are proportional to A 2 , B 2 and AB respec- 
tively. 

Explicitly, they are given by 

(r!|C(0,k 1 )C(0,k 2 )C(0,k 3 )C(0,k 4 )|r!) SB A2 = (2^) 3 <5( 3 )(K)ll? 1 (fc 1 2/fc)P c 2 (fc)P c (fc 1 2) + 2permutations, (C2) 



2 Strictly speaking, Eq. (135 I t is only valid when all momentum vectors have similar magnitude as only in this case we can treat, for 
instance, H as a constant. However the momentum dependence of (135 1 is such that T^^* a ' goes to zero in the squeezed limit and this 
is consistent with (I B 1 ft 
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(r!|C(0,k 1 )C(0,k 2 )C(0,k 3 )C(0,k 4 )|r!) S£ Ai 5 = (2ir) 3 6 i3) (K)\v 2 (k 12 /k)V^k)V c (k 12 ) + 2 permutations, (C3) 

(fi|C(0,k 1 )C(0,k 2 )C(0,k 3 )C(0,k4)|O) SB B2 = {2TT) 3 S ( - 3) (K)^V 3 (k 12 /k)V^{k)V c (k 12 ) + 2 permutations, (C4) 

where we have defined three T>i(x) functions, with i — 1, 2, 3, as 

9a; 4 (3a; 5 + 36a; 4 + 176a; 3 + 432a; 2 + 528a; + 512) 



X>i(a:) 
V 2 (x) 



2 7 (a; + 2) 6 

-3a; 4 (13a; 7 + 156a: 6 + 748a; 5 + 1696a; 4 + 1072a; 3 - 4160a; 2 - 11712a; - 8192) 

2 7 (a; + 2) 6 : 



x 4 

^3 (x) = ^ X (103a: 9 + 1236a; 8 + 5276a; 7 + 5632a; 6 - 27052a; 5 - 95248a; 4 - 69376a; 3 

v ' 2 9 (a; + 2) 6 v 

+ 149056x 2 + 281792a; + 131072), (C5) 



where k\ 2 = ky/2(l + cos 6*3), fci3 = ky/2(l + cos 2 ) and kn — k^j2(l + cos6*i). k denotes the common amplitude 
of all four momentum vectors kj, with j = 1, ... ,4 and the rescaled power spectrum is V^(k) = H 2 / (Aec s k 3 ) that 
should be evaluated at horizon crossing. 



http:/ /map. gsfc.nasa.gov/ 

WMAP, E. Komatsu et al, Astrophys. J. Suppl. 180, 330 (2009), 0803.0547. 
K. M. Smith, L. Senatore, and M. Zaldarriaga, (2009), 0901.2572. 
J. M. Maldacena, JHEP 05, 013 (2003), |astro-p h/02 10603] 
|http: / / www.rssd.esa.int / index. php?project=Pla nck| 

A. P. S. Yadav and B. D. Wandelt, Phys. Rev. Lett. 100, 18130 1 (2008), 0712.114 8. 

A. D. Linde and V. F. Mukhanov, Phys. Rev. D56, 535 (1997), |astro-ph/9610219| 

N. Bartolo, S. Matarrese, and A. Riotto, Phys. Rev. D65, 103505 (2002), h ep-ph/0112261 
F. Bernardeau and J.-P. Uzan, Phys. Rev. D66, 103506 (2002), hep-ph/0207295 

F. Bernardeau and J.-P. Uzan, Phys. Rev. D67, 121301 (2003), astro-ph/0209330] 

G. Dvali, A. Gruzinov, and M. Zaldarri aga, Phys. Rev. D 69, 023505 (2004)7 |astoo-ph/0303591~j 
P. Creminelli, JCAP 0310, 003 (2003), [astFo-ph/0306122| 

M. Alishahiha, E. Silverstein, and D. Tong, P hys. Rev. D 70, 123505 (2004) , |hep-th /0 404084 
A. Gruzinov, Phys. Rev. D71, 027301 (2005), [astro-ph/0406129| 

K. Enqvist, A. J okinen, A. Mazumdar, T. Multamaki, and A. Vaihkonen, Phys. Rev. Lett. 94, 161301 (2005), 
|astro-ph/0411394| 

A. Jokinen and A. Mazumdar, JCAP 0604, 003 (2006), [astro-ph/0512 368 

D. H. Lyth, JCAP 0511, 006 (2005), |astro-ph/0510443| 

M. P. Salem, Phys. Rev. D72, 123516 (2005), |astro-ph/0511146| 

D. Seery and J. E. Lidsey, JCAP 0701, 008 (2007), |astro-ph/0611034| 

M. Sasaki, J. Valiviita, and D. Wands, Phys. Rev. D7 4, 103003 (2006), [a"stro-ph/0607627| 

K. A. Malik and D. H. Lyth, JCAP 0609, 008 (2006) , |astr"o-ph/0604387| 

N. Barnaby and J. M. Cline, Phys. Rev. D73, 10 6012 (2 006), astro-ph/060 1481| 

L. Alabidi and D. Lyth, JCAP 0608, 006 (2006), |asteo-ph/0604569j 

X. Chen, M.-x. Huang, S. Kachru, and G. Shiu, JCAP 0701, 002 (2007), hep-th/0605045 
X. Chen, M.-x. Huang, and G. Shiu, Phys. Rev. D74, 121301 (2006), hep-th/0610235 
X. Chen, R. Easther, and E. A. Lim, JCAP 0706, 023 (2007), astro-ph/0611645 
L. Alabidi, JCAP 0610, 015 (2006), [astFo-ph/06046lT] 

C. T. Byrnes, M. Sasaki, and D. Wands, Phys. Rev. D74, 123519 (2006), [astFo-ph/0611075| 
T. Suyama and M. Yamaguchi, Phys. Rev. D77, 023505 (2008), 0709.2545. 

F. Arroja and K. Koyama, Phys. Rev. D77, 083517 (2008), 0802.1167. 

F. Arroja, S. Mizuno, and K. Koyama, JCAP 0808, 015 (2008), 0806.0619. 

D. Langlois, S. Renaux-Petel, D. A. Steer, and T. Tanaka, Phys. Rev. Lett. 101, 061301 (2008), 0804.3139. 
D. Langlois, S. Renaux-Petel, D. A. Steer, and T. Tanaka, Phys. Rev. D78, 063523 (2008), 0806.0336. 

M. Sasaki, Prog. Theor. Phys. 120, 159 (2008), 0805.0974. 

C. T. Byrnes, K.-Y. Choi, and L. M. H. Hall, JCAP 0810, 008 (2008), 0807.1101. 
C. T. Byrnes, K.-Y. Choi, and L. M. H. Hall, JCAP 0902, 017 (2009), 0812.0807. 



10 



[37] B. Dutta, L. Leblond, and J. Kumar, Phys. Rev. D78, 083522 (2008), 0805.1229. 
[38] A. Naruko and M. Sasaki, Prog. Theor. Phys. 121, 193 (2009), 0807.0180. 
[39] T. Suyama and F. Takahashi, JCAP 0809, 007 (2008), 0804.0425. 
[40] X. Gao, JCAP 0806, 029 (2008), 0804.1055. 

[41] H. R. S. Cogollo, Y. Rodriguez, and C. A. Valenzuela-Toledo, JCAP 0808, 029 (2008), 0806.1546. 
[42] Y. Rodriguez and C. A. Valenzuela-Toledo, (2008), 0811.4092. 

[43] K. Ichikawa, T. Suyama, T. Takahashi, and M. Yamaguchi, Phys. Rev. D78, 023513 (2008), 0802.4138. 
[44] C. T. Byrnes, JCAP 0901, 011 (2009), 0810.3913. 

[45] S. Li, Y.-F. Cai, and Y.-S. Piao, Phys. Lett. B671, 423 (2009), 0806.2363. 

[46] D. Langlois, F. Vernizzi, and D. Wands, JCAP 0812, 004 (2008), 0809.4646. 

[47] C. Hikage, K. Koyama, T. Matsubara, T. Takahashi, and M. Yamaguchi, (2008), 0812.3500. 

[48] M. Kawasaki, K. Nakayama, T. Sekiguchi, T. Suyama, and F. Takahashi, JCAP 0811, 019 (2008), 0808.0009. 

[49] X. Gao and B. Hu, (2009), 0903.1920. 

[50] Y.-F. Cai and H.-Y. Xia, (2009), 0904.0062. 

[51] D. Langlois, S. Renaux-Petel, and D. A. Steer, (2009), 0902.2941. 

[52] X. Gao, (2009), 0904.4187. 

[53] Q.-G. Huang, (2009), 0903.1542. 

[54] E. Silverstein and D. Tong, Phys. Rev. D70, 103505 (2004), [hep^th/03102211 

[55] S. Kecskemeti, J. Maiden, G. Shiu, and B. Underwood, JH EP 09, 076 (2006) |hep-th/0605189] 
[56] J. E. Lidsey and D. Seery, Phys. Rev. D75, 043505 (2007), |astro-ph/ 0610398 

[57] D. Baumann and L. McAllister, Phys. Rev. D75, 123508 (2007), hep-th/061 0285] 

[58] R. Bean, S. E. Shandera, S. H. Henry Tye, and J. Xu, JCAP 0705, 004 (2007), |hep-th/0702107| 
[59] J. E. Lidsey and I. Huston, JCAP 0707, 002 (2007), 0705.0240. 

[60] H. V. Peiris, D. Baumann, B. Friedman, and A. Cooray, Phys. Rev. D76, 103517 (2007), 0706.1240. 
[61] T. Kobayashi, S. Mukohyama, and S. Kinoshita, JCAP 0801, 028 (2008), 0708.4285. 
[62] L. Lorenz, J. Martin, and C. Ringeval, JCAP 0804, 001 (2008), 0709.3758. 
[63] S. Bird, H. V. Peiris, and D. Baumann, (2009), 0905.2412. 



[64] L. Boubekeur and D. H. Lyth, Phys. Rev. D73, 021301 (2006), astro-ph/0504046 

[65] L. Alabidi and D. H. Lyth, JCAP 0605, 016 (2006) , |astro-ph/05 10441 1 

[66] N. Kogo and E. Komatsu, Phys. Rev. D73, 083007 (2006), astro-ph/0602099 

[67] D. Babich, P. Creminelli, and M. Zaldarriaga, JCAP 0408, 009 (2004), astro-ph/0405356 

D. Seery, J. E. Lidsey, and M. S. Sloth, JCAP 0701, 027 (2007), [aitrb-ph/0610210j 

[69] D. Seery, M. S. Sloth, and F. Vernizzi, JCAP 0903, 018 (2009), 0811.3934. 

[70] J. Garriga and V. F. Mukhanov, Phys. Lett. B458, 219 (1999), [hep-th/9904176] 

[71] C. Armendariz-Picon, T. Damour, and V. F. Mukhanov, Phys. Lett. B458, 209 (1999), |hep-th/9904075l 

[72] A. J. Christopherson and K. A. Malik, (2008), 0809.3518. 

[73] S. Weinberg, Phys. Rev. D72, 043514 (2005) , |hep^th/050 6236 

[74] P. Adshead, R. Easther, and E. A. Lim, (2009), 0904.4207. 

[75] P. Creminelli and M. Zaldarriaga, JCAP 0410, 006 (2004), [astFo-ph/0407059] 

[76] C. Cheung, A. L. Fitzpatrick, J. Kaplan, and L. Senatore, JCAP 0802, 021 (2008), 0709.0295. 

[77] X. Chen, B. Hu, M.-x. Huang, G. Shiu, and Y. Wang, (2009), 0905.3494. 



